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In coalescing neutron star binaries, r-modes in one of the stars can be resonantly excited by the 
gravitomagnetic tidal field of its companion. This post-Newtonian gravitomagnetic driving of these 
modes dominates over the Newtonian tidal driving previously computed by Ho and Lai. To leading 
order in the tidal expansion parameter R/r (where R is the radius of the neutron star and r is the 
orbital separation), only the I = 2, |m| = 1 and |m| — 2 r-modes are excited. The tidal work done 
on the star through this driving has an effect on the evolution of the inspiral and on the phasing of 
the emitted gravitational wave signal. For a neutron star of mass M, radius R, spin frequency /spin, 
modeled as a F = 2 polytrope, with a companion also of mass M, the gravitational wave phase 
shift for the m = 2 mode is ~ 0.1 radians (7?/lOkm)''(A'//1.4M0)-i°/^(/spin/lOO Uzf/^ for optimal 
spin orientation. For canonical neutron star parameters this phase shift will likely not be detectable 
by gravitational wave detectors such as LIGO, but if the neutron star radius is larger it may be 
detectable if the signal-to-noise ratio is moderately large. The energy transfer is large enough to 
drive the mode into the nonlinear regime if /spin ^ 100 Hz. For neutron star - black hole binaries, 
the effect is smaller; the phase shift scales as companion mass to the -4/3 power for large companion 
masses. The net energy transfer from the orbit into the star is negative corresponding to a slowing 
down of the inspiral. This occurs because the interaction reduces the spin of the star, and occurs 
only for modes which satisfy the Chandrasekhar-Friedman-Schutz instability criterion. 

A large portion of the paper is devoted to developing a general formalism to treat mode driving in 
rotating stars to post-Newtonian order, which may be useful for other applications. We also correct 
some conceptual errors in the literature on the use of energy conservation to deduce the effect of 
the mode driving on the gravitational wave signal. 



I. INTRODUCTION AND SUMMARY 
A. Background and motivation 

One of the most promising sources for ground-based 
gravitational wave observatories such as LIGO [l| and 
VIRGO 3 are coalescences of compact binary systems 
[sj. The first searches for these systems have already 
been completed [3| , and future searches will be more sen- 
sitive. For neutron star-neutron star (NS-NS) binaries, 
the most recent estimate of the Galactic merger rate is 
83lg6^ Myr"^ [1, Q; extrapolating this to the distant 
Universe using the method of Ref. [tI yields a rate of 
28l22yr"^ within a distance of 200 Mpc. The range of 
next-generation interferometers in LIGO is expected to 
be ~ 300 Mpc [|. Coalescing binaries may well be the 
first detected sources of gravitational waves. 

The gravitational waveforms will carry a variety of dif- 
ferent types of information [1, Q: First, in the early, low 
frequency phase of the signal, the binary can be treated 
to a good approximation as consisting of two spinning 
point masses. The corresponding waveforms give detailed 
information about the deviations of general relativity 
from Newtonian gravity, and they have been computed to 
post-3. 5-Newtonian order 9]. Second, toward the end of 
the inspiral (frequencies of order several hundred Hertz, 
see Fig. [Ij the internal degrees of freedom of the bodies 
start to appreciably influence the signal, and there is the 
potential to infer information about the nuclear equation 
of state. This has prompted many numerical simulations 



of the hydrodynamics of NS-NS mergers using various 
approximations to general relativity (see, eg, Ref. [Toj 
and references therein). Equation of state information 
can also potentially be extracted from the waves energy 
spectrum and from the NS tidal disruption signal 
for neutron star-black hole (NS-BH) binaries p^ . 

A third type of information potentially carried by the 
gravitational waves is the effect of the internal structure 
of the bodies on the early, low frequency (10 Hz < / < 
100 Hz) portion of the signal, via the excitation of in- 
ternal modes of oscillation of one of the neutron stars 
by the tidal gravitational field of its companion. In this 
regime the waveform's phase evolution is dominated by 
the point-mass dynamics, and the perturbations due to 
the internal modes are small corrections. Nevertheless, 
if the accumulated phase shift due to the perturbations 
becomes of order unity or larger, it could impede the 
matched-filtering-based detection of NS-NS signals [§]. 
Alternatively the detection of a phase perturbation could 
give information about the neutron star structure. For 
these reasons the excitation of neutron star modes has 
been studied in detail. 

A given mode can be treated as a damped, driven har- 
monic oscillator, described by the equation 

q + jq + <g - A{t) cos[m(j){t)]. (1.1) 

Here g is a generalized coordinate describing the mode, 
cjin is the mode frequency, 7 is a damping constant, (j){t) 
is the orbital phase of the binary, m is the azimuthal 
quantum number of the mode, and A(t) a slowly varying 
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FIG. 1: Noise curves /irms(/) = sj fSh{f) for initial and next- 
generation LIGO interferometers [l^l (shown in thin lines) 
together with the signal strength hc{f) for an inspiralling bi- 
nary of two IAMq neutron stars at a distance of 200 Mpc 
(thick line) [13 ]. The signal will terminate near the frequency 
of the innermost stable circular orbit, which is ~ 850 Hz 
as shown here if the neutron stars have radii of 10 km and 
are modeled as F = 2 polytropes 15]. Resonantly excited 
modes can give rise to small but potentially measurable cor- 
rections to the phase of the waveform in the early portion 
10 Hz < / < 100 Hz of the signal. 



amplitude. As the inspiral of the binary proceeds, the 
orbital frequency lo = (j) gradually increases with time. 
Sufficiently early in the inspiral, therefore, the mode will 
be driven below resonance: muj <C uj^^. Because the 
inspiral is adiabatic we also then have A/ A <^ uj^ and 
cj/cj <C Using these approximations Eq. (Il.ip can be 
solved analytically in this pre-resonance regime, and one 
finds that the total energy absorbed by the mode up to 
time t is 

ft 



, uj^^/m^ + ijj(t 



This energy absorption causes the orbit to inspiral 
slightly faster and changes the phase of the gravitational 
wave signal. 

Mode excitations can have three different types of ef- 
fects on the gravitational wave signal: 

• Dissipative: There is a phase perturbation due to 
the second term in Eq. (|1.2p . corresponding to en- 
ergy dissipated by the mode. This is a cumulative 
effect. It is dominated by the I = 2,m = ±2 funda- 
mental modes of the neutron star, and its magni- 
tude depends on the size of the damping constant 7, 
which is determined by the fluid viscosity. This ef- 
fect was examined by Bildsten and Cutler [l6| , who 
showed that the phase perturbation is negligible for 
all physically reasonable values of the viscosity. 

• Adiabatic: There is a phase perturbation due to 
the first term in Eq. (|1.2p . This corresponds to 



the instantaneous energy present in the mode, as 
it adjusts adiabatically to sit at the minimum of 
its slowly-evolving potential. The perturbation to 
the waveform is again dominated by the I = 2,m — 
±2 fundamental modes Ij/Tj. This effect has been 
studied analytically in Refs. fT6l.jT8l . IT9l . [20l. [2l| and 
numerically in Refs. [TtI. [2^. [2311241 [25|. It can also 
be inferred from sequences of 3D numericalquasi- 
equilibrium models of NS-NS binaries [TTI. [26l. [27j. 
The phase perturbation grows like uj{t)^^^ and is of 
order a few cycles by the end of inspiral. It may be 
marginally detectable for some nearby events [2^ . 

• Resonant: The frequencies of internal modes of a 
neutron star are typically of order Wm ~ \/M/R^, 
where M and R are the stellar mass and radius. 
(We use units with G — c = 1 throughout this pa- 
per). For a NS-NS binary, this is of order the orbital 
frequency at the end of the inspiral, ~ 1000 Hz. 
Therefore, most modes are driven at frequencies be- 
low their natural frequencies throughout the inspi- 
ral, and are never resonant. However, some mode 
do have frequencies with tJni ^ \/ M / R^^ and 
these modes are resonantly excited as the gradu- 
ally increasing driving frequency muj sweeps past 
LOm- During resonance the effect of damping is neg- 
ligible, and the width Aw of the resonance is de- 
termined by the inspiral timescale. The energy de- 
posited in the mode is thus enhanced compared to 
the equilibrium energy [the first term in Eq. (jl.2p ] 
by the ratio of the inspiral timescale to the orbital 
period. This can be a significant enhancement. 

Resonant driving of modes in NS-NS binaries has been 
studied in Refs. [B, HI, [H, HI, [H, H^, [pi, and in the 
context of other binary systems in Refs. |32l . [3^ . [s^, [3^. 
One class of modes with suitably small frequencies are 
(/-modes; however the overlap integrals for these modes 
are so small that the gravitational phase perturbation is 
small compared to unity [T3, [SO] . Another class are 
the / and p-modes of rapidly rotating stars, where the 
inertial-frame frequency Wm can be much smaller than the 
corotating-frame frequency ^ yjM/R-^. Ho and Lai jsTj 
showed that the phase shifts due to these modes could be 
large compared to unity. However, the required NS spin 
frequencies are several hundred Hz, which is thought to 
be unlikely in inspiralling NS-NS binaries. 

A third class of modes are Rossby modes (r-modes) 
[3^ . 113] , for which the restoring force is dominated by the 
Coriolis force. For these modes the mode frequency 
is of order the spin frequency of the star, and thus can be 
suitably smaU [10 Hz < uj ,„/(27r) < 100 Hz]. Ho and Lai 
[3TI ] computed the Newtonian driving of these modes, and 
showed that the phase shift is small compared to unity. 

In this paper we compute the phase shift due to the 
post-Newtonian, gravitomagnetic resonant driving of the 
r-modes. Normally, post-Newtonian effects are much 
smaller than Newtonian ones. However, in this case 
the post-Newtonian effect is larger. The reason is that 
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the Newtonian driving acts via a coupling to the mass 
quadrupolc moment 



Qij = J (fxp{xiXj - -^(JtjX^) 



(1.3) 



of the mode, where p is the mass density. For r modes, 
this mass quadrupole moment vanishes to zeroth order 
in the angular velocity of the star; it is suppressed by 
a factor ~ ^B?' jM. By contrast, the gravitomagnetic 
driving acts via a coupling to the current quadrupole 
moment 

Sij = ^ J d^^P [(x X + (x X \-)jXi\ , (1.4) 

where v is the velocity of the fluid. This current 
quadrupole moment is nonvanishing to zeroth order in 
O, and the gravitomagnetic coupling is therefore en- 
hanced over the Newtonian tidal coupling by a factor 
of ^ VP'B? /M . This enhancement factor can exceed the 
usual factor ~ v'^ /c^ by which post-Newtonian effects are 
suppressed, where w is a velocity scale. It is for this same 
reason that the radiation-reaction induced instability of 
r-modes in newly born neutron stars [ssl . [s^ is domi- 
nated by the gravitomagnetic coupling and not the mass 
quadrupole coupling. 



B. Effect of resonance on gravitational wave signal 

We now discuss the observational signature of a mode 
resonance in the gravitational wave signal. We denote 
by <^(t) the phase of the waveform, which is twice the 
orbital phase (f){t). We denote by $pp(t) the phase that 
one obtains from a point particle model, neglecting the 
effects of the internal modes of the stars. To the leading 
post-Newtonian order, this point-particle phase is given 
from the quadrupole formula by the differential equation 

m 
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(1.5) 



where w = is the orbital angular velocity, and Mt and 
fj, are the total mass and reduced mass of the binary. 

Consider now the phase <I>(i) including the effect of 
the mode. We neglect the phase shift caused by the adi- 
abatic response of the mode at early times before the 
resonance [the first term in Eq. (II. 2p ] , as this phase shift 
is small compared to the effect of the resonance. We 
also neglect the cumulative phase shift due to damping. 
Therefore at early times we have $(t) = $pp(i). The 
duration of the resonance is short compared to the inspi- 
ral timescale (see Sec. II CI below). After the resonance, 
the mode again has a negligible influence on the orbit of 
the binary, and Eq. p.Sp applies once more. However, 
the two constants of integration that arise when solving 
Eq. (jl.Sp need not match between the pre-resonance and 



post-resonance waveforms. Therefore, the phase can be 



written as 



Kn(t 



ppW 

ppV^ -r M) 



^0 ^ ^ir( 



(1.6) 



Here is the time at which resonance occurs, and tics 
is the duration of the resonance.^ The effect of the res- 
onance is thus to cause an overall phase shift A<I>, and 
also to cause a time shift At in the signal. If energy is 
absorbed by the mode, then the inspiral proceeds more 
quickly and At is positive. 

It turns out that the two parameters At and A$ char- 
acterizing the effect of the resonance are not independent. 
To a good approximation they are related by 



A$ = $pp(to)At, 



(1.7) 



as argued by Reisenegger and Goldreich [29|, and as de- 
rived in detail in Sec. IIVI below. Inserting this relation 
into Eq. (|1.6p and expanding to linear order in the small 
parameter At yields for the late-time phase 



m = *pp(<) + 'i'ppW - ^pp(to) 



At. 



(1.8) 



The physical meaning of the condition (|1.7p is thus 
that the resonance can be idealized as an instantaneous 
change in frequency at t = to with no corresponding in- 
stantaneous change in phase. Rewriting the phase ((TT 
in terms of A<I> using (|1.7p and using <I>pp — 2uj gives 

"^(t) 



- 1 



OJQ 



A$, 



(1.9) 



where luq — uitQ) is the orbital frequency at resonance. 
The phase perturbation due to the resonance therefore 
grows with time after the resonance, and can become 
much larger than A<I> when w 3> wq- 

One can alternatively think of the phase correction as 
being non-zero before the resonance, and zero afterwards 
[2^. This is what would be perceived if one matched 
a template with the portion of the observed waveform 
after the resonance. In other words, if we define a point- 
particle waveform phase $pp(t) with different conven- 
tions for the starting time and starting phase via 



$pp(t) = $pp(t + At)-A$, 
then we obtain 



(1.10) 



-pp 
?pp 



u(t) 



1 



A$ 



to <C — ti-c 
to 3> trns- 



(1.11) 



^ The sign of the A$ term in Eq. III.6I I is chosen for later conve- 
nience. 

^ At intermediate times \t~to \ ~ trcs, the phase 3>(t) is not related 
in a simple way to <t>pp{t). However, since the duration of the 
resonance is typically a few tens of cycles (short compared to 
the ~ 10* cycles of inspiral in the detector's frequency band), 
the phase perturbation in this intermediate stage will not be 
observable. 
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The phase correction <J> — <I>pp is now a hnear function 
of frequency, interpolating between A<i> at a; = to zero 
at resonance lo — uq. Since this phase correction never 
exceeds A$, we see that the detectabihty criterion for 
the resonance effect is A$ > 1. Thus, the phase shifts 
^ A$ predicted by Eq. ()1.9p are fictitious in the sense 
that they are not measurable. 



C. Order of magnitude estimates 



(jl.lSp is small compared to unity, so that the force on 
the mode is in phase with the mode's natural oscillation. 
This gives ties ^ V^orb^rr, and we get 



A<I> aiji 



(1.16) 



Finally, we use the scaling t„ ^ r^/{iJLM^) for the ra- 
diation reaction time given by the quadrupole formula, 
where ji is the reduced mass. This gives 



We now turn to an order of magnitude estimate of the 
mode excitation and the phase shift A$, both for the 
Newtonian driving studied previously and for the gravit- 
omagnetic driving. 

We consider a star of mass M and radius R with a 
companion of mass M' . The orbital angular velocity at 
resonance is wq! this is also the mode frequency (up to 
a factor of the azimuthal quantum number m which we 
neglect here). The separation r of the two stars at reso- 
nance is then given by cjq ~ Mt /r^ , where Mt = M + M' 
is the total mass. 

Let the magnitude of the tidal acceleration be atid, and 
let the duration of the resonance be tres- The mode can 
be treated as a harmonic oscillator of mass ~ M and 
frequency ajQ. During the resonance, the mode absorbs 
energy at the same rate as a free particle would, so the 
total energy absorbed is 



(1.12) 



Now the gravitational wave luminosity Eg^ is to a good 
approximation unaffected by the resonance. Therefore 
the loss of energy AE from the orbit decreases the time 
taken to inspiral by an amount At = AE/Egy^. From 



Eq. (|1.7p the corresponding phase shift is A(f> ^ At /to 
where torb ~ l/i^o is the orbital period. This gives 



A$ 



t,,AE 



(1.13) 



where i?orb is the orbital energy and t^.^ is the radiation 
reaction timescale. This phase shift is just the energy 
absorbed divided by the energy radiated per orbit. 

Next, using E^rh ~ MM' /r and the formula pTT2)) for 
AE gives 



A$ ^ -a. 



t^ t 



tid" 



to 



(1.14) 



The resonance time ij-es is the geometric mean of the or- 
bital and radiation reaction times torb and t-^^. This fol- 
lows from the fact that the orbital phase near resonance 
can be expanded as 

(^^^0 + Looit -to) + ^{t- tof + ..., (1.15) 

from the definition of the radiation reaction timescale t^^. 
The mode is resonant when the quadratic term in Eq. 



A$ 



M'^i'^M^ 



(1.17) 



We now consider various different cases. For the New- 
tonian driving of a mode via a tidal field of multipole of 
order we have 



Otid 



M' 



i-i 



(1.18) 



[The leading order, quadrupolar driving is I = 2.] Insert- 
ing this into Eq. (|1.17p and eliminating r using 



,2^3 



Mt gives 



A$ 



M3 



Mj \M/R^ 



{l-5)/3 



(1.19) 



which agrees with previous analyses [29|, |30| for g and / 
modes. 

For the Newtonian driving of r-modes, the leading or- 
der driving is / = 3 rather than I = 2 [3i| . Also the tidal 
acceleration atid should be multiplied by the suppression 
factor ojqR'^/M discussed in Sec. II Al above, and the final 
phase shift p.lOp should be multiplied by the square of 
this factor. This gives in the equal mass case M = M' 



10, ,10/3 



A$ 



rt LO, 



M20/3 



(1.20) 



Now the frequency mtJo of the dominant, I ^ m = 2 r- 
mode is related to the spin frequency /spin of the neutron 
star by 2ujq = 87r/spijj/3, so we obtain 



A$ 



7->10 j-10/3 
_ J spin 

' M20/3 



(1.21) 



where /c is a dimensionless constant of order unity. Equa- 
tion (|1.2ip agrees with the results of Ho and Lai [31| . 
who show that k « 0.4 for a F = 2 polytrope, and for the 
I = m = 2 r-mode.'^ This gives 



A$r- 



modc.N 



= 8 X 10-^i?i«/:{'/3'M-f /^ (1.22) 



^ See equation (4.13) of Ref. [sil . specialized to the equal mass 
case and maximized over the angle /3 between the spin of the 
star and the orbital angular momentum. 
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where the subscript N denotes Newtonian and we have 
defined R = 10i?iokm, M = IAM^Mq, and /spi„ = 
100/,iooHz. 

For the gravitomagnetic driving of r-modes, the tidal 
acceleration is given by 



atid 



M'R 



^spin^oib ■ 



(1.23) 



Here the first factor in brackets is tidal acceleration (|1.18p 
for Newtonian quadrupolar / — 2 driving. Since the in- 
teraction is gravitomagnetic, it is suppressed by two pow- 
ers of velocity relative to the Newtonian interaction. One 
power of velocity is associated with the source of the grav- 
itomagnetic field, and is the orbital velocity Uorb ~ '^o''- 
The second power of velocity is that coming from the 
Lorentz-type v x B force law, and is the internal fluid 
velocity in the neutron star due to its spin, Uspin ^ ujqR. 
Combining these estimates with Eq. p.l7p and eliminat- 
ing luq in favor of /spi„ gives 



n>4 j"2/3 
J spin 

M'M'^mI'^ ' 



(1.24) 



In the equal mass case M = M' this gives 



A$r-modc,PN ~ 3.4 -R^o/'ioo^iTi"^'- (1-25) 
Here the subscript PN denotes post-Newtonian. 
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FIG. 2: The magnitude of the gravitational-wave phase shift 
A"l> due to the r-mode resonance as a function of the spin 
frequency /spin of the neutron star, for both the Newtonian 
and gravitomagnetic driving, for a binary with masses M = 
M' — IAMq, and for a neutron star radius of i? = 10 km. 
The star is modeled as a F = 2 polytrope. The mode plotted 
is i = 2, m = 2 in each case. The angle tp between the spin 
and the orbital angular momentum has been chosen in each 
case to maximize the phase shift; tl) is 34.4° in the Newtonian 
case and 60° in the post-Newtonian case. The resonance oc- 
curs at an orbital frequency of 4/spin/9 in the Newtonian case 
and at 4/spin/3 in the post-Newtonian case. The correspond- 
ing gravitational wave frequencies are 8/spin/9 and 8/spin/3 
respectively. It can be seen that the gravitomagnetic driving 
dominates. 



D. Results and implications 

The fact that the order-of-magnitude estimate p.25p 
of the phase shift is of order unity, and thus potentially 
detectable in the gravitational wave signal via matched 
filtering is the primary motivation for the computa- 
tion of this paper. We compute the dimcnsionless coef- 
ficient in Eq. p.24p . and find for the case of a F = 2 
polytrope for the I = 2, m = 2 mode [cf. Eq. (|5.42p be- 
low, multiplied by two to convert from orbital phase shift 
A(j) to gravitational wave phase shift A$] 

A$r-tnodc.PN - -0.06 Rjo^M^l^/^ . (1.26) 

This is somewhat smaller than the estimate (ll.25p . Fig- 
ure O compares this result to the Newtonian phase shift 
(fr22)) . In TableHwe show the phase shift A$22 = 2A022 
that will occur at coalescence for the five known double 
neutron star binaries.^ For B1913-I-16 we used the known 
value tp = 22° of the angle between the pulsar spin and 



The orbital motion of these five systems is eccentric today, but 
radiation reaction will have circularized the orbits by the time 
resonance occurs. It is thus consistent to use our circular-orbit 
results for these binaries. 



System 


spin period 


M 


M' 


phase lag A $2 2 


B1913-fl6 


59.0 


1.441 


1.387 


-0.005 


B1534-I-12 


37.9 


1.339 


1.339 


-0.009 


B2127+11C 


30.5 


1.358 


1.354 


-0.010 


J0737-3039 


22.7 


1.34 


1.25 


-0.013 


J1756-2251 


28.5 


1.40 


1.18 


-0.011 



TABLE I: List of gravitational wave phase lag parameters 
that will occur for the Z = 2, m = 2 resonance for the 
known double neutron star binaries when they coalesce. For 
B1913-I-16, the known value for the misalignment angle tp 
(22°) was used; for the others, the average value of the an- 
gular factors was used. The pulsar mass AI and companion 
mass M' are in solar masses and the pulsar spin period is in 
milliseconds. The parameter values were taken from Ref. • 

the orbital angular momentum [4]| . and for the remain- 
ing binaries for which ip is unknown we took an average 
of the angular factors appearing in the phase lag (j5.37bp . 
We assumed a F = 2 polytrope equation of state and a 
radius of 10 km for all pulsars. 

The phase shift (|1.26p is sufficiently small that it is un- 
likely to be detectable in the gravitational wave signal for 
most detected binary inspirals. Even for the most favor- 
able conceivable neutron star parameters of M = 1.0 



6 



and R = 15 km, the phase shift is ~ 0.8 radians. An 
analysis of the detectabihty of mode resonances indicates 
that the minimum detectable value of A$ is A<i> ~ 2 for 
/spin = 100 Hz at a signal to noise ratio of S/N — 10 
[43|. Therefore the minimum signal to noise ratio nec- 
essary for detection of the effect is between ^ 30 and 
300, depending on the neutron star parameters. It is 
conceivable that the effect could be detected in some of 
the closest detected binary inspirals with advanced inter- 
ferometers. 

One of the interesting features of the interaction is that 
the phase shift is negative, corresponding to energy be- 
ing transferred from the star to the orbit rather than the 
other way around. This arises because the star is rotating 
and the interaction reduces the spin of the star, tapping 
into its rotational kinetic energy.^ In Sec. IIVB] below we 
show that this sign of energy transfer is only possible 
for modes which satisfy the Chandrasckhar-Friedman- 
Schutz mode instability criterion [44, 45] of being retro- 
grade in the corotating frame and prograde in the inertial 
frame. 

Although the r-mode resonance may not be detectable, 
it is nevertheless probably the strongest of all the mode 
resonances in the early phase of the gravitational wave 
signal. The energy transfer corresponding to the phase 
shift (fL^ is from Eq. ([13^ 

AE ^ (1.4 X 104« erg) i?^o/sioo, (1-27) 

which is fairly large. The corresponding dimensionless 
amplitude a of the r-mode in the notation of Rcf. [46] is 

a^O.SRljsWoM^-l (1.28) 

This amplitude is sufficiently large that nonlinear mode- 
mode interactions should be important, as in the case 
of unstable r-modes in newly born neutron stars and in 
low mass X-ray binaries ^46,] . In those contexts the mode- 
mode coupling causes the mode growth to saturate at a '--^ 

Here however the r-mode is 
driven on much shorter timescales, and the mode-mode 
coupling likely does not have time to operate efficiently. 
It is possible that th e type of nonlinear effects studied 
numerically in Refs. [53, [5J] (generation of differential 
rotation, shocks, wave breaking) could occur. 



E. Organization of this paper 

The remainder of this paper is organized as follows. In 
Sees . nil IhII and HVl we analyze resonant driving of modes 
in Newtonian gravity, building on the previous studies 



of Refs. m, [13, mi, H Hi. The time delay parameter 
At can be derived using energy conservation which we 
discuss in Appendix [b1 but the relation p.7p between 
the phase shift parameter A$ and At requires an anal- 
ysis of the orbital equations of motion, which we give in 
Sec. IIVI In Sec. [Vl and Appendix [XI we show that the 
gravitomagnetic driving of modes can be computed to 
leading order using Newtonian stellar perturbation the- 
ory supplemented by a gravitomagnetic external acceler- 
ation, and we deduce the phase shifts generated by the 
driving of r-modes. 



II. NEWTONIAN MODE-ORBIT COUPLING 

In this section we describe a model dynamical system 
that describes the Newtonian excitation of modes in in- 
spiralling binary systems. We generalize previous treat- 
ments [sj, [35I ] to allow for stellar rotation. The dynami- 
cal system will actually be general enough to encompass 
the post-Newtonian driving of r-modes, as we discuss in 
Sec. [Vl below. 



A. Mode expansions in rotating Newtonian stars 

We start by reviewing the evolution equations for mode 
amplitudes in rotating Newtonian stars, following the 
treatment in Ref. [13]. We assume that the background 
star is uniformly rotating with angular velocity CI and 
mass M. We denote by {x,t) the inertial frame coor- 
dinates and by {x,t) the co-rotating frame coordinates. 
These coordinate systems are oriented so that $1 is in 
the z direction. Below we will use spherical polar coor- 
dinates (r, 9s,4>s) for X] the corresponding inertial frame 
coordinates are (r,9s,(f)s) where 

<j)s=4's+^t. (2.1) 

A mode of a rotating star is a pair {^{x),uj) of a mode 
function ^(x) and a rotating-frame frequency w, such 
that the rotating-frame Lagrangian displacement given 
by 

^(S,i)=|(S)e— * (2.2) 

is a solution to the linearized hydrodynamic equations. 
We label the different modes with w < ^ by an index j; 
the jth mode is {^j{x),ujj). It can then be shown that 
the Lagrangian displacement ^{x,t) and its time deriva- 
tive ^{x,t) can be expanded in the phase-space mode 



^ This follows from the fact that at fixed angular momentum, the 
energy of a star is minimized at uniform rotation. Therefore the 
only way the total energy of a rotating star can be reduced is via 
a reduction of the spin. 



In a rotating star, to obtain a complete set of modes for the phase 
space mode expansion, one can either choose all the modes with 
oj > or all the modes with uj < 0. Here we choose the latter 
(so the relevant r-modes have m > 0); the opposite convention 
was used in Ref. |47|| . 
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expansion [47 



c*it) 



(2.3) 



In rotating stars, one cannot in general find a basis of 
modes which are real, so the expansion coefficients Cj{t) 
are complex in general. 

If there is an external acceleration acy±{x^t) acting on 
the star, then the equation of motion for the jth mode is 

M 



where the inner product is defined by 



(2.5) 



and p{x) is the background density. The constant hj is 
given by 

h,^{S,,,2ujjij+2iVLxi,). (2.6) 

We normalize the modes using the convention 

(C„^,)=Afi^^ (2.7) 

where R is the stellar radius. 

An equivalent second-order form of the equation of mo- 
tion can be obtained by differentiating Eq. (|2.4[) and elim- 
inating Cj in the result using Eq. (j2.4p . This gives 

+w|cj = fj - iujjfj, 

where fj{t) = acxt)/&j- This is the form of the equa- 
tion of motion discussed in the Introduction, cf. Eq. (jl.ip . 



B. Mode-orbit coupling equations 

We now consider a binary system of a star of mass M, 
whose modes are excited, and a companion of mass AI' 
which we treat as a point mass. We denote the total 
mass M + M' by Mt, and the reduced mass MM' /Mt 
by ji. We denote by a; = (r, 6*, (/>) the displacement of 
M' with respect to M in polar coordinates. We orient 
these polar coordinates so that the motion is confined to 
the equatorial plane 6 — 7r/2. The equations of motion 
describing this system are 



rcff 



modes I diss 



2f </> = a" 



odes 



(2.8a) 
(2.8b) 
(2.8c) 



Here Cj and ujj are the amplitude and rotating-frame fre- 
quency of the jth mode. The right hand sides of Eqs. 



(I2.8bp and (|2.8cp are the components of the total relative 
acceleration a"i°dcs _|_ ^diss ^.j^g orthornomal basis e^;, 
e?, where a"i°'i'=>5 describes the back-reaction of excited 

modes onto orbital motion and a^^^^ is the back-reaction 
of emission of gravitational radiation onto orbital mo- 
tion. The latter acceleration term is obtained from the 
Burke- Thorne radiation-reaction acceleration 15511 



^diss 



(i,x) 



5 dt^ 



(2.9) 



where Qij is the quadrupole moment of the binary and 
TF means trace-free part. In evaluating this expres- 



(2.4) sion we include only the contribution Qi 



flXiXj to 



the quadrupole from the orbital degrees of freedom, ne- 
glecting the contribution from the modes. To evaluate 
the time derivatives of Qij we use the Newtonian or- 
bital equations of motion with the mode coupling terms 
omitted. We also neglect the effect of gravitational wave 
damping on the modes themselves. This yields the com- 
ponents of the dissipative acceleration 



and 



l6^iMt 



SpMt 



2 J,2 



9r^ - 6r 



2 i2 



4Mt 
3r 



2Mt 



(2.10) 



(2.11) 



The function fj in Eq. (|2.8a|) describes the coupling 
between mode j and the tidal gravitational field of the 
companion object. For Newtonian mode-orbit coupling, 
it can be derived as follows [Slj. Insert into the right 
hand side of Eq. (|2.4p the acceleration a^xt — — V$oxt, 
where the external potential is 



$ext = -V 



M' 



X — x(t)\ 

oo /' 
l'=2 m' = ~l' 



I -im' 4>{t) 



it) 



i+i 



(2.12) 



Here the coefRcient Wi'm' is given by Eq. (2.2) of Ref. 
[sH . and x{t) = (r(t), 7r/2, 0(i)) is the location of the 
companion. The spherical harmonics that appear in 
Eq. (|2.12p are related to to spherical harmonics of the 
inertial- frame coordinates {r,6s,(t>s) which are aligned 
with the spin axis of the star by 



Y.'^il'Yi'AOs^^s), (2.13) 



where V^i^, is the Wigner P-function [31|. Finally use 



mm' 

the relation (12. ip between the inertial and corotating 
frame coordinates 6^ and 6^. The result is 



fj — Qj,l'm' 

I 



(t)+i 



ifit 



ity 



+1 



(2.14) 
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where rrij is the azumithal quantum mimber of the mode 



^j, and 



iM' . 



i j,l'm' 



-Wvrn'-D. 



(2.15) 

Note that there are two different azimuthal quantum 
numbers in the final expression (j2.14p . the quantum num- 
ber rrij of the mode which controls the number of factors 
of the phase e'^* , and the quantum number m' which con- 
trols the number of factors of the phase e"*"^^*-*. These 
two quantum numbers coincide when the spin axis of the 
star is aligned with the orbital angular momentum, but 
not in general. 

For the Newtonian case, in the absence of dissipa- 
tion, the dynamical system (j2.8p can be derived from 
the Hamiltonian 



2 pI fiMt 



2fi ' 2^r2 

(2.16) 

where pr = fir, p^ — fJ.r'^cf), together with the Poisson 
brackets {cj^Ck} = 0, {cj, c^} — —iSjk/bj. This allows 
us to read off the orbital acceleration a'^o^cs ^j^g ^j^g 
modes: 



^modcs 



odes 



fir 



3 flA. ^3 fIJ. 



(2.17b) 



From the formula (|2.14p we see that a given mode j is 
driven by a sum of terms labeled by {l',m'), the effects 
of which superpose linearly. For the remainder of this 
section and in Sees. IIIII and IIVI we will focus attention 
on a particular mode j and on a particular driving term 
{l',m'). Using the notations Aj{r) = \Qj,i'm'\/r'' '^^ and 
e'^J — Qj.i'm' /IQjj'ni'l, the coupling function (|2.14p and 
the orbital acceleration p.l7bp can be written as 



1] 

^modcs 



fj — [7.(i)] e-™'</'(t)+i™3f^t+iX3 ^ 



(2.18a) 



2m'b,A,[rit)] 
fir{t) 



—im' <p(t)-\-imj Q,t-\-ixj 



(2.18b) 



Although these formulae (|2.18p were derived for Newto- 
nian tidal driving, we shall see in Sec. fVl below that they 
also apply to gravitomagnetic driving in the vicinity of 
a resonance, with appropriate choices of the functions 
Aj{r) and phases Xj- 

C. New orbital variables 

In order to analyze the effects of the mode coupling 
on the orbital motion it is useful to make a change of 
variables, from r{t), (j){t), r{t), (j){t) to new variables p(t) , 



(f>{t), e{t) and (/>o(i) that characterize the Newtonian el- 
lipse that is instantaneously tangent to the orbit. Here 
p is the instantaneous semilatus rectum and e is the in- 
stantaneous eccentricity. These variables are defined by 
the equations 



P 



(2.19a) 
(2.19b) 
(2.19c) 



where = cos(0 — (po) and — sin(0 — (j)Q). In terms 
of these variables the orbital equations become 




2p- 



,3/2 



^(1 + ec^) ^ 
/lit, 



p- 



,3/2 



■(l + ec^)^ 



(2.20a) 
(2.20b) 



2VM(1 + ec^) 

X (3e -I- 4c0 -I- ec2,/,)a ; 



(2s^ -I- es20)af 



(2.20c) 



-VP 



2eVMt(l + ec0) 

X (e -I- 2c0 -I- ec2^)af 



(4s0 -I- es24>)a^ 



(2.20d) 

Here = cos[2(0 - 0o)], = sin[2((/) - (/jq)], and 
af = af°'^'''' + ap'' and = aj°'^'''' + a'^^'''^ are the compo- 
nents of the total acceleration in the radial and tangential 
directions. 



D. The circular approximation 

In this paper we restrict attention to situations where 
the eccentricity e is negligibly small before, during and 
after the mode resonance. We will justify this assump- 
tion in Sec. IIV CI below by estimating the magnitude of 
the eccentricity generated during the resonance. Using 
this approximation we can simplify the set of equations 
as follows. We can drop Eqs. (|2.20c[) and (j2.20dp , and 
we simplify the remaining equations (|2.8ap . (|2.20ap . and 
(j2.20bp using e = 0. Using the expressions (|2.10p and 
(|2.1ip for the dissipative components of the acceleration. 
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and the mode-orbit coupling terms (I2.18p . now gives 

Cj+iuojc = ^^-(p)e-™>(*)+'"^-"*+*x^ (2.21a) 



Am'hjAj{p)^ 



,3/2 



x3fi 
0. 



^ — im <f}{t)+imjQt~\-iXj 



F 



, (2.21b) 

(2.21c) 



Note that only the (j) component of the acceleration con- 
tributes in this approximation. 



E. The no-backreaction approximation 

When the amplitude of the mode coupling Aj is small 
we can solve the system (|2.2ip of equations by neglecting 
the backreaction of the mode excitation on the orbital 
motion. More precisely, we (i) solve the orbital equa- 
tions (j2.21bp and (j2.21cp for the zeroth order solutions 
p(t) and 0(i), neglecting the mode coupling; (ii) insert 
those zeroth order solutions into Eq. (|2.21ap to compute 
the evolution of the mode amplitude c(i); and (iii) insert 
that mode amplitude c{t) into Eqs. (|2.21bp and (|2.21cp to 
obtain the linearized perturbations Sp{t) and S(p{t) to the 
orbital motion. Steps (i) and (ii) are carried out in Sec. 
mil where we extend previous results in the vicinity of 
resonance by Rathore, Blandford and Broderick 3^ 3^ . 
Step (iii) is carried out in Sec. IIVI We justify the use of 
this no-backreaction approximation in Sec. IIVDI 



III. MODE AMPLITUDE EVOLUTION 

A. Zeroth order solutions 

We start by reviewing the zeroth order solutions of the 
orbital evolution equations l|2.21bp and (|2.21cp that apply 
in the limit Aj ^ of no modal coupling. We denote by 
to the time when mode j enters resonance, and by loq 
the orbital angular velocity at resonance. We define the 
radiation reaction timescale at resonance Ur to be 



We define the dimensionless time parameters 



A = 



(3.3) 



(time from resonance in units of the resonance time) , and 

t-to 



tr 



(3.4) 



(time from resonance in units of the radiation reaction 
time) . These variables will be useful in our computations 
below. We also define the dimensionless small parameter 



^rcs 

0.015 



32 /6 



3 V 5 
M 



5/6 



1.2 Af, 



/o 



10 Hz 



5/6 



(3.5) 
(3.6) 



We will use e as an expansion parameter throughout our 
computations below. It follows from these definitions 
that e is also the ratio of the orbital and resonance times, 
up to a constant factor 



3e' 



We also have 



eA. 



(3.7) 



(3.8) 



Using these notations, the zeroth order solutions can 
be written as 



m = 0o + ^[l-(l-r)^/' 
Lo{t) = = c.o(l-r)-3/8, 
Pit) - po(l - r)i/^ 



(3.9a) 

(3.9b) 
(3.9c) 



Here 0o, ujq and po are the values of (f), lu and p at reso- 
nance, related by 



Wo 



3/2 

Po 



(3.10) 



tr 



8/3 



157 s 



M 



1.2 M, 



-5/3 



/o 
10 Hz 



-8/3 



(3.1) 



where A4 — fi'^^^M^^^ is the chirp mass and /o = 
u!o I (27r) . We also define the resonance timescale fjcs via 
[cf. Sec. HO] above] 



(96Al5/3^^^/V5 



-1/2 



= 2.5s 



M 



1.2 M, 



-5/6 



/o 



10 Hz 



-11/6 



(3.2) 



B. Mode ampHtude evolution 



In the rest of this section we will solve Eq. (|2.21ap for 
the mode amplitude analytically using matched asymp- 
totic expansions. We will divide the inspiral into three 
regimes: an "early time" regime before resonance, the 
resonance regime, and a "late time" regime after reso- 
nance. We will compute separate solutions in these three 
regimes and match them in their common domains of 
validity. 

The evolution equation (|2.21ap for the mode amplitude 
can be written as 

Cj + iujjCj — Aj{p) exp [imjnt — im'(l){t) + ixj] , (3.11) 
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where and (j){t) are given by Eqs. (|3.9p . It can be seen 
from this equation that resonance occurs at an orbital 
angular frequency luq given by 



m'ujQ = mjfl + LOj. 



(3.12) 



Note that the right hand side of this equation is just the 
inertial-frame frequency of the mode. We will restrict 
attention to cases for which the sign of m' is the same 
as the sign of rajfl + Uj, so that the solution c^o of Eq. 
p.l2p is positive and resonance does occur. 



using Eq. (|3.7p and the resonance condition (|3.12p gives 



cKA) = 



■ exp[i(itr — o^jtrosA — to'A^/2)] 



1 + 



A2 



(3.15) 



where Ur = Xj ~ ™'0o + Wj^lio- Here the error term 
0(1/ A^) comes from the error term in Eq. (|3.13p . while 
the second term O(eA^) arises from expanding the expo- 
nential in the expression p.l4ap .^ The expression (I3.15P 
is therefore vaHd for 1 < |A| < e^^/^. 



1. Early time solution 

Well before the resonance the amplitude and frequency 
of the forcing term on the right hand side of Eq. p. lip are 
changing on the radiation reaction timescale iir , which is 



much larger than the mode period 



and the pe- 



riod 



of the tidal forcing. Therefore to a first 



approximation we can neglect the time-dependence of 
Aj{t) = Aj[p{t)] and Lj{t) = <j>{t). This gives the fol- 
lowing approximate solution to Eq. (j3.1ip : 



i[mjQ — m'uj{t) + LOj 
X J 1 + 



(3.13) 



Here the superscript e denotes "early time", and we have 
used the initial condition that the mode excitation van- 
ishes as i ^ — oo. The error estimate in Eq. (I3.13P can be 
obtained by writing the mode amplitude as [t) + 5d^ (t) 
with d'it) given by Eq. (P+^ . substituting into Eq. (jXllD 
and solving for the linearized correction Sd^^t) using the 
same method of neglecting the time dependence of uj{t) 
and A{t). The error terms show that Eq. p.lSp is valid 
in the regime — cjol S> From the explicit solu- 

tion (j3.9bp for ijj{t) this condition can also be written as 
|r| > e, or, using Eq. ([3^ . as |A| > 1. (Recall that A 
and T are negative before resonance). 

For later use it will be convenient to specialize the ex- 
pression (|3.13p for the early time solution to the regime 
|i — to| <C irr or equivalently |A| ^ 1/e. Expanding the 
zeroth order expressions (|3.9p about t — to, writing the 
result in terms of the rescaled time A and using Eq. p.7p 
gives 



8A 1 

4>o + — + -X' + 0{eX'), (3.14a) 

cjo [l + 3eA/8 + 0(e2A2)] , (3.14b) 

Aoj[l + 0(eA)]. (3.14c) 



2. Late time solution 

The late-time solution after resonance can be found 
using the same method as for the early time solution. 
The only difference is that we need to add an arbitrary 
solution of the homogeneous version of Eq. (|3.1ip . which 
will be determined by matching onto the solutions in the 
resonance and early time regions. The late-time solution 
can be written as 



j\^^.(j-^^[i™,jnt-im' (j>(t)+ixj\ 

i[mj^ — m'uj{t) + LUj] 



X U + o 



(w - UJoY 



(3.16) 



where Bi and xi ^^"^ constants parameterizing the homo- 
geneous solution and the superscript "1" denotes "late 
time". As before, this solution is valid in the regime 
|i — tol ^ ^res or |A| ^ 1; in the late time regime A is pos- 
itive. Also as above we can further expand the solution 
in the near-resonance regime 1 <C |A| <C e^^^'^ as 



4(A) 



3m'wo eA 



exp[«(ur — Wj^rcsA — m'A /2)] 



0(eA3 



(3.17) 



3. Resonance time solution 

Finally we turn to the resonance regime |t — to| ~ ^res 
or |A| ~ 1. In this regime we can expand the differen- 
tial equation (|3.1ip using the expansions (|3.14p and the 



where Aqj is the value of Aj at resonance. Su bstitu ting f^xs^ the amplitude is negli gible in comparison i 

these expansions into the early time solution (|3.13P and error term from the expansion I l3.14al l of the phase 



Since A S> 1 here the O(eA) error term from the expansion 
l|3.14c| l of the amplitude is negligible in comparison to the O(eA^) 
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formulae (13.31) and (13.71) to obtain 



which can be expressed in terms of the Fresnel integrals 



+ iujjt.csc'j = Aojtrcs X [1 + O(eA)] 



7 "^rcs , ^,..vj (3.18) 

^gi[«r-iJjtroBA-m'A^/2+0(eA^)] 



f C(z) 1 




cos 








I J 





du, 



(3.19) 



The solution to this equation consists of a homogeneous 

solution which will be determined by matching to the as discussed by Rathore, Blandford and Broderick [3^ . 



early-time solution, together with a particular solution 



J 



cKA) 



1 . 




The solution is 



r 



(3.20) 



Here the superscript "r" denotes "resonance" and Br and 
Xr are constants parameterizing the homogeneous part of 
the solution. The error term O(eA^) in Eq. p.20p comes 
from the correction term in the argument of the exponen- 
tial in Eq. (|3.18p . Note also that the error terms gener- 
ated from the errors in the amplitude expansion p.l4cp 
scale as eexp(zA^) or e^A and can be neglected. It fol- 
lows from Eq. (|3.20p that this resonance solution is valid 
in the regime |A| iT^I"^ . 

For later use it will be convenient to specialize the ex- 
pression p.20p for the resonance time solution to the 
regime |A| 3> 1. In this regime the arguments of the 
Fresnel integrals are large, and we can use the asymp- 
totic formulae 



C(z) 1 ^ sgn(z) \_ 

S{Z) J 2 TTZ 




This gives 



{^z'/2) + 0{z-^). 

(3.21) 



rZZsgn(A)e*("'--'^^*'-^-^s"('"')'^/2) 
V 2|r7i'| 



_|__^gi(nr-WjtrcsA-m'A^/2) 

m'A 

[1 + 0(A-3) + 0{tX^)] + Bre"'^'^^'''^^+^-\ 

(3.22) 

This form of the resonance solution is valid in the regime 

1 « lAI « £-1/2. 



4-. Matching computation 



We now match the expressions p.lSp for the early time 
solution and (|3.22p for the resonance solution within their 
common domain of validity 1 <C |A| <C e~^/^ with A < 0. 
For both expressions the dominant error terms in Cj(A), 
which are obtained by multiplying the fractional errors 



by the corresponding expressions for Cj(A), are of order 
0{\^) and 0(e-iA-3). Using Eq. (p!7)l one can show 
that expression p.lSp matches the second term inside 
the square brackets in Eq. (|3.22p . Demanding that the 
remaining terms in Eq. (j3.22p cancel one another deter- 
mines the constants Br and Xr'- 



Br 



2 TO' 



— sgn(m')7r/2) 



(3.23) 

Here we have used the fact that sgn(A) = —1 before 
resonance. In a similar way we match the expressions 
p.l7p for the late time solution and (|3.22p for the reso- 
nance solution within their common domain of validity 
1 < |A| < £-1/'"^ with A > 0. This yields Bi = 2Br and 
Xi = Xr, or 



Bie'^' 



271 



-Ao,,i„se^("'-^sn(m')-/2) + o{e)] 



(3.24) 

The size of Bi implies that at late times when |w — Wj| > 
LUj , the homogeneous term in the late time solution p.l6p 
is larger than the particular solution by a factor 1/e. 
Therefore in this regime the late-time solution is a freely 
oscillating mode. 



IV. PERTURBATION OF ORBITAL MOTION 

In this section we solve the orbital evolution equations 
(j2.21bl - I2.21cp , treating the tidal terms as linear per- 
turbations on top of the quasi-circular binary inspiral 
motion. We look for solutions of the form 



Pit) 

m 



Spit), 



(4.1a) 
(4.1b) 



where p^'^^ (t) and (j)^^^ (t) are the zeroth order inspiral so- 
lutions (|3.9cp and p.9ap . It will be sufhcient to specialize 
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to the regime |r| <C 1 near resonance. Linearizing Eqs. 
(|2.21bl - I2.21cp yields the following evolution equations 
for the perturbations 5p and 5(j): 



Sp - 



4irr(l-r) 



6p = s{t), 
Sp = 0. 



—irn' 4){t)-\-irnj Q.t-\-ixj 



2p^/2(l-r)V8 
Here the source term s{t) is 

s{t) = 2so^oj5R 
where 

2to'6jVpo 

So = 7=- 



A. Evolution of semi-latus rectum 



(4.2a) 
(4.2b) 



[l + 0(eA)], 
(4.3) 



(4.4) 



Equation (I4.2a|) is a first order differential equation and 
its solution is easily found to be 



Sp^ I sit') f 1 '^^'^ 



(4.5) 



where r' = t'/^rr- In evaluating integral (|4.5p . we use the 
early time approximation (|3.13p for the mode amplitude 
c{t) for ~oo < t < to — tni, the resonance time approx- 
imation (|3.20p for to — tn^ < t < to + and the late 
time approximation (|3.16p for t > io + ^m- We denote 
by s'^{t), s'^{t), and s'(t) the corresponding approxima- 
tions to the driving term ()4.3p . and by Sp'^{t), 5p'^{t) and 
5p^{t) the corresponding approximate solutions for 5p in 
the three regimes. We choose the parameter tm govern- 
ing the times at which we switch from one approximation 
to the next to be 



at„e 



4/5 



(4.6) 



where cr is a dimensionless constant of order unity. The 
corresponding values oi t — [t — io)/irr and A = r/e are 
Tm = cre^/^ and Am = cre^^/'"'. The choice (|4.6I) of scaling 
with e of maximizes the overall accuracy, as can be 
seen from the scalings of the error estimates in the mode 
amplitude solutions fXTS)) . (jXTT)) and ((2221) ■ The final 
results are independent of the choice of a. 



Performing integral ()4.5p is straightforward, albeit a 
little tedious. The result of interest to us here is the 
expansion of 6p^{t) in the regime 1 <C A <C e^^^^. This is 



V(i) = ^Mi^f- 



2 

res 



i-i ^ 

AV Trim' 



|m'|A2 



TT 
4 



27r|m'|A2 



(4.7) 



B. Phase evolution 

We obtain the phase perturbation Scj) by inserting into 
Eq. (j4.2bp the perturbation Sp of the semi-latus rectum 
and integrating. We first discuss the leading order phase 
shift and then discuss the magnitude of the corrections. 

To leading order we take 5p{t) = for t < to, and for 
t > to we use the expression for Sp{t) given by the first 
term in the brackets in Eq. (|4.7p . This approximate form 
of Sp{t) satisfies for t > to the homogeneous version of the 
equations of motion (|4.2p . It coincides at |t| ^ 1 with the 
homogeneous solution obtained by taking to to — At in 
the zeroth-ordcr solution given by Eqs. (|3.4p and (|3.9cp 
and linearizing in At: 



Sp{t) = - 



PoAt 



4(1 - Tf/H„ 



(4.8) 



By expanding Eq. (|4.8p in the regime 1 ^ A <C e~^^^, 
and comparing with the first term in Eq. (|4.7p . we can 
read off the parameter At: 



At- 



Svr 



|m'|po 



so^^,,4.in- [1 + 0(e)] . (4.9) 



We discuss below the 0(e) error term. 

The corresponding solution for the phase perturbation 
5(j} is also a solution of the homogeneous system of equa- 
tions, and therefore must be of the form discussed in the 
introduction, parameterized by the time delay parameter 
At and an overall phase shift A0: 

0(i) = (/)(°) (t) + 5(i,{t) = {t + At) - Acj). (4.10) 



In this approximation the relation (|1.7p between A0 and 
At follows from the continuity of S<l){t) at t = to, which 
follows from Eq. (|4?2b)) : 



A(f) = ujoAt. 



(4.11) 



Combining this with the expression (|4.9p for At and the 
formula (|4.4p for sq gives the phase shift formula^ 

A0 = 167r sgn(m') ' '"'^ " [1 + 0(e)] . (4.12) 



* The gravitational wave phase shift A<I> used in the introduction 
is related to the orbital phase shift A<f> used here by A$ = 2A0. 
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We next discuss the sign of the phase shift A0. We 
shall see later [cf. Eqs. l|iTT|) and ([BO]) ] that the sign 
of the phase shift is the same of the sign of the net en- 
ergy transferred from the orbit to the star. The sign of 
the parameter bj in Eq. (|4.12[) coincides with the sign of 
LUj for incrtial modes, and probably also in general [sgI ]. 
Therefore we obtain 

sgn(A0) — sgn{m'LL)j) — sgn{mjn + LUj) sgn^LUj) 

— sgn{Q + ujj/mj) sgn{Lj J /nij), (4.13) 

where on the first line we have used the resonance condi- 
tion (|3.12p . Since O > by convention it follows that 
A(f> is negative if and only if both Uj/mj is negative 
(the mode is retrograde in the corotating frame) and 
Uj I ruj + Q is positive (the mode is prograde in the in- 
ertial frame), i.e., the Chandrasekhar-Friedman-Schutz 
mode instability criterion [H, is satisfied. 

Consider now the magnitude of the corrections to the 
leading order result (|4.12p . First, there is a nonzero 
phase perturbation in the early time regime, generated by 
the perturbation 6p. This phase perturbation grows like 
1/|A|, but is smaller than the phase shift ()4.12p by a fac- 
tor of ~ e when the growth saturates at |A| ~ 1. Second, 
there is a phase perturbation which is generated by in- 
serting the second and third terms in the approximation 
for Sp'^{t), and integrating using Eq. ([42bl) . This 
phase perturbation is also is smaller than Acj) by a factor 
of ~ e. Third, there are fractional corrections of order 
0(e) due to higher order contributions to the matching 
parameters Br, Xr, Bi and xi indicated by the error terms 
in Eqs. (|3.23p and (|3.24p . Finally, there are phase pertur- 
bations generated at late times by the coupling between 
the freely oscillating mode and orbital motion. These can 
also be shown to be smaller than the leading order A(f> 
by a factor of at most ~ e. 



C. Validity of circular approximation 



When solving the orbital evolution equations (|2.20ap 
and Ij2.20bp we set the eccentricity e to zero. We now 
determine the domain of validity of this approximation. 

First, in the absence of mode coupling, there is a non- 
zero instantaneous eccentricity e ~ ^orb/^rr ~ due to 
the gradual inspiral. To compute this eccentricity we 
start with the orbital equations of motion (j2.8bp and 
p.8cp with the mode coupling terms dropped, and solve 
using a two timescale expansion. The result is that the 
zeroth order inspiral solution (j3.9p is accurate up to frac- 
tional corrections of order O(e^): 

rit) = ro(l-T)i/4[l + 0(e4)], (4.14a) 

(t,{t) = / dtLUoil - t)--^/« [1 + Oie^)] . (4.14b) 



Pit) = ro(l-r)i/4[i + o(,4)] ^ 
eit) ^ 



using the definition (|2.19p gives 

(4.15a) 

-(l-r)-'V8 [1 + 0(^2)], (4.15b) 
m = j dtuj^{l - r)-3/8 [1 + 0(^4)] ^ (4 15^) 
Mt) = W) - 3^2] [1 + 0(e')] . (4.15d) 

We now denote this solution by p(")(<), e(")(i), 0(°)(t) 

and (p^Q^ (t), and linearize the equations (|2.20p to solve for 
the perturbations Sp{t), 6e{t), 54>{t) and 5(f)Q{t) generated 
by the mode coupling. We retain only terms that are 
zeroth order in e. The resulting equation for 5e{t) is 



5e 



Mt 



mode 



(4.16) 



ode 



coupling, given by Eqs. (|2.17ap and (|2.18ap . We now use 
the expression (|3.20p for the mode amplitude in the res- 
onant regime, simplify using the expansions (|3.14p . and 
in the resulting expression drop all the terms that do not 
accumulate secularly. The result is 



Seit) = 6e^,^J{^\rn^\X), 
where the function J'(A) is given by 



(4.17) 



^(A) = / dXi - sin(A72) 



-cos(AV2) 



S 



A 



-'4 



This function satisfies J'(A) ^ as A ^ 00, from the fact 
that the Fresnel functions C and S are odd and using the 
approximate formulae (|3.2ip . Therefore in this approxi- 
mation the final eccentricity generated by the resonance 
vanishes, in agreement with a previous result of Rathore 
[Eq. (6.104) of Ref. ^ speciahzed to e = 0]. 

The quantity i5eniax is, up to a factor of order unity, the 
maximum eccentricity achieved during the resonance. It 
scales as 



e'Acj). 



(4.18) 



Transforming to the variables p{t), e{t), (j){t) and 4>o{t) 



This eccentricity is transient and lasts only for a time 
^ tics- From the equations of motion (|2.20aP and (j2.20bp 
we can estimate the corrections to the parameters A(j) and 
At characterizing the resonance that are generated by 
this transient eccentricity. We find that the corrections 
to both A(j) and LUoAt are of order eA(j) [(5emax times the 
number of cycles ~ 1/e of resonance] or smaller, and thus 
can be neglected. 

We note that this computation is valid only in the 
regime Acj) <C 1 for which Se <^ e^^^ throughout the res- 
onance. In the regime Acj) ^ 1, both e{t) and cjio{t) are 
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significantly perturbed away from their zeroth order so- 
lutions, and it is not a good approximation to evaluate 
the right hand side of Eq. I|4.16p using the zeroth order 
solutions as done here. We conclude that the circular 
approximation is valid for this paper, since At/) < 1 for 
all the cases we consider. 



D. Validity of no-backreaction approximation 

The no-backreaction approximation is valid when the 
phase perturbation S(j){t) accumulated over the resonance 
is small compared to the phase accumulated over the 
same period in the absence of mode couphng [s^l, i.e., 
when 



T 



adapted to the star whose modes are being driven. We 
can specialize this coordinate system so that (i) the post- 
1-Newtonian mass dipole of the star vanishes, so that 
the origin of coordinates coincides with the star's cen- 
ter of mass; (ii) the angular velocity of the coordinate 
system as measured locally using Coriolis-type accelera- 
tions vanishes (the angular velocity with respect to dis- 
tant stars will not vanish due to dragging of inertial 
frames); and (iii) the normalization of the time coor- 
dinate is such that the I = piece of the Newtonian 
potential outside the star is of the form 4> = —M/r with- 
out any additional additive constant (sgI . [57| . In such 
body-adapted reference frames (also called local asymp- 
totic rest frames) the effect of the external gravitational 
field on the internal dynamics of the star can be param- 
eterized, in post-l-Newtonian gravity, by a set of grav- 
<C 1. (4.19) itoelectric tidal tensors £i-^i2...iiit) and gravitomagnetic 



This ensures the smallness of the backreaction fractional 
corrections to the mode driving terms on the right hand 
side of Eq. ([XTT]) . From Eq. ((j!2b)) we have 



ty. 



5/2 

Po 



Sp{to). (4.20) 



Using Eq. (|4.5p to estimate Sp{to) now gives 



(4.21) 



Comparing this with the formula (I4.12p for the resonance 
phase shift and using Eq. (|3.7p gives 

T - e^A(f>. (4.22) 



The estimate (|4.22l) shows that in general there is a 
non-empty regime where the no-backreaction approxima- 
tion is valid, T <C 1, and where in addition the phase shift 
due to the resonance is large, A(/) ^ 1. For the r-modes 
studied here, by combining the numerical estimate (|1.26p 
of A^ with the numerical estimate p.Sp of e wc obtain 



T . 



10-3 Af-f^^^/^ 



ft /■J n4 
JslOO^lO- 



(4.23) 



Thus the no-backreaction approximation is valid for the 
mode driving considered here. 



V. APPLICATION TO GRAVITOMAGNETIC 
TIDAL DRIVING OF ROSSBY MODES 

A. Overview 



tidal tensors B. 



(t), for Z > 2. These tensors are 



symmetric and tracefree on all pairs of indices, and are 
invariant under the remaining post-l-Newtonian gauge 
freedom [5^. 

In Newtonian gravity the gravitoelectric tidal fields are 
just the gradients of the external potential <I>cxt evaluated 
at the spatial origin: 



1 



{I - 2)1 dx'^ . . . dx^^ 



iO,t). (5.1) 



In post-l-Newtonian gravity the definition of the tidal 
tensors is more complicated, since the field equations are 
nonlinear and so the potentials can not be written as the 
sum of interior and exterior pieces as in the Newtonian 
case. The post-l-Newtonian definitions are discussed in 
detail in Refs. [H, \5^.^ The tidal tensors Si^i 



and 

can also be defined in the more general context 



of a small object in an arbitrary background metric [58l |. 

In this section we will compute the driving of the r- 
modes by the leading order gravitomagnetic tidal tensor 
Bij{t). The effects of the higher order tensors Bi^,,,ii{t) 
for I > 3 are suppressed by one or more powers of the 
tidal expansion parameter i?/r, where R is the radius of 
the star and r is the orbital separation. Therefore we will 
neglect these higher order tensors. 

A key point about the gravitomagnetic tidal tensors is 
that they vanish identically at Newtonian order. There- 
fore when computing the driving of the r-modes, it is suf- 
ficient to use Newtonian-order stellar perturbation theory 
supplemented by post-l-Newtonian tidal driving terms. 
This is sufficient to compute the mode driving to post- 
l-Newtonian accuracy, and constitutes a significant sim- 
plification.^*^ 



In this section we compute the parameters At and A<I> 
that characterize the perturbation (|1.6p to the gravita- 
tional wave phase due to the resonant post-l-Newtonian 
tidal driving of r-modes. 

In our computations we will use a harmonic, con- 
formally Cartesian post-l-Newtonian coordinate system 



We follow here the notational conventions of Ref. [58l for the 
tidal tensors fij^.i, (t) and Bj, (t). A different notational con- 
vention is used in Refs. [s^. |57(| . where these tensors are denoted 
-Gii...i,{t)/(«-2)! and -3Hi-^,„i^ {t)/[2{l+l){l-2)\] respectively. 
Note that the validity of this approximation depends on the 
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For the Newtonian mode driving discussed in the pre- 
vious sections, we used the Hamiltonian (|2.16p to deduce 
driving terms p.l7p in the orbital equations of motion. 
Here, for the post-Newtonian case, we will instead com- 
pute those driving terms directly, by evaluating the cur- 
rent quadrupole Sij induced by the r-mode and using the 
post-l-Newtonian equations of motion including multi- 
pole couplings derived in Ref. [s^. We will verify that 
these driving terms take the same form (j2.18bp as in the 
Newtonian case and as analyzed in Sees. Illll and IIVI 

The remainder of this section is organized as follows. 
We review the properties of r-modes in Sec. IV Bl Section 
IV CI and Appendix [X] compute the mode-orbit coupling 
term responsible for driving the mode, as well as the 
resonant response of the modes. Next in Sec. IV Dl we 
analyze the effect of the mode on the orbit and derive 
the parameters At and A$, using the results of Sec. IIVI 



B. Rossby modes 

In this section we review the relevant properties of 
Rossby modes (r-modes). We will use the notations of 
Sec. Ill Al above, in particular the inertial frame coordi- 
nates are x — {r,0s,(j>s) and the coordinates that co- 
rotate with the star are x = (r, 9s,4>s) with (f>s = 4>s -\-^t. 
In the remainder of this section we will drop the mode in- 
dex i for simplicity, and we will write the rotating-frame 
mode frequency ujj as a;i„. 

For r-modes, the rotating frame mode frequencies are 

M 

2mVl . , 

where 17 = |r2| and m are the usual spherical harmonic 
indices. The mode eigenfunctions are [3g| 

£,ir.{x) = -^^M= s X VY,Ue, 4>) [1 + 0{vu')] , 
V H' + 1) 

(5.3) 

where is a spherical harmonic and fim{r) is a real 
radial mode function. Here w = ri/i7cr is the angular 
velocity in units of the break-up angular velocity ficr = 
a/ M/R^. The mode normalization condition (|2.7p can 
be written in terms of the radial mode function /;,„ (r) as 

drr^p{r)fira{rf ^MR". (5.4) 



Using the mode function (|5.3p in the definition (|2.6p of 
the constant b yields 

h = LUn.MR'^. (5.5) 



choice of coordinate system. It is valid for the coordinate sys- 
tems used here, but would not be valid for a coordinate system 
{t, x') obtained from the standard global inertial-frame harmonic 
coordinate system (t, i') by a Newtonian- type transformation of 
the form t = t, = — z*{i). 



We define a dimensionless coupling parameter /;.,„ for the 
mode by 



Iim — 



1 



r^firnir)pir)dr; (5.6) 



this parameter will arise in our computations below. 

We next express the mode function (|5.3|1 in inertial 
coordinates {x,t), in which it simply acquires a factor of 



■^—iraVtt 

ii.m{x) = j==fi^ X X VYi^{0, ^) [1 + Oiw'')] . 

V H' + 1) 

(5.7) 

For each m we define the symmetric tracefree tensors 
X™2...s, by the equation [H 

YUO,cf)^y::s....s,n''n^' ...n^\ (5.8) 

where n* = /r is the unit radial vector. 



C. Computation of force terms and resonant 
response of modes 

We now turn to the analysis of the mode amplitude 
driving terms (^/m, ctoxt)- We restrict attention to modes 
^;„i with m > 0; the to < driving terms are just the 
complex conjugates of the m > terms. In order to 
derive the form of the external acceleration aoxt(ic,i) to 
be added to the Newtonian perturbation equations, we 
consider, temporarily, post-l-Newtonian stellar pertur- 
bation theory. The argument which follows is a slightly 
more rigorous version of the argument given in Appendix 
A of Ref. [131. 

In post-l-Newtonian gravity in conformally Cartesian 
gauges, the metric is expanded in the form 



2e^^ + 2e^($2 



i,) + 0{e'')]{dt/e) 



-[2e^Ci+0{e^)]dx\dt/e) 



0{e^)]dx'dx^ . 



(5.9) 



Here $, ■0 and C, are the Newtonian potential, the post- 
Newtonian scalar potential and the gravitomagnetic po- 
tential, respectively. The quantity £ is a formal expan- 
sion parameter which can be set to unity at the end of 
the calculation; it can be thought of as the reciprocal of 
the speed of light. 

We consider a uniformly rotating background star, 
characterized by a pressure poj density po, fluid 3- velocity 
Vq, and by gravitational potentials $0i V'Oi and Co- We 
assume that the star is subject to an external perturbing 
gravitomagnetic tidal field Bij{t). The Hnear response 
of the star can be parameterized in terms of the Eule- 
rian perturbations Sp, Sp, Sv, (5$, Sip, and S(^, which sat- 
isfy the linearized post-l-Newtonian hydrodynamic and 
Einstein equations. The solution is determined by the 
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boundary conditions on the gravitational potential per- 
turbations at large r, which are 



6^Ij{x, t) 



(5.10a) 

(5.10b) 
(5.10c) 



as r — > cx). [See, for example, Eqs. (3.5) of Ref. [56| |. 
where we drop all the tidal tensors except for Hij = 
—2Bij.] Next, the hnearized equation satisfied by SC in 
harmonic gauges is 



SC - IGttpoSv' - 16Trv'„Sp = 0. 



(5.11) 



We can neglect the second and third terms in this equa- 
tion for the reason explained in Sec. IV Al above: the per- 
turbations dp and Sv will be proportional to Bij and thus 
of post-l-Newtonian order, so the corrections to SC they 
generate will be of post-2-Newtonian order. Solving Eq. 
(15. lip without the second and third terms and subject to 
the boundary condition (|5.10ap yields 

5Q{x,t) = ^e,jkBki{t)x'xK (5.12) 

Also, since the perturbation has no Newtonian part, we 
have 



= 0. 



(5.13) 



Next, we take the field equation for ip, the post- 
l-Newtonian continuity equation, and the post-l- 
Newtonian Euler equation in harmonic gauge. We lin- 
earize these equations about the background solution 
(PO)Po, 'I'Oj ^0: V'Oj Co), use Eq. (|5.13p and drop all terms 
of the type discussed above that give rise to post-2- 
Newtonian corrections. The results are 



and 



6p + V ■ {SpvQ + pqSv) = 0, 
V5p 



(5.14) 



(5.15) 



Po 



5v — — (fo • ^)5v — {5v ■ V)t;o 
where 

Qoxt = -5C + f X (V X 50 



(5.16) 



Finally, we can replace the post-l-Newtonian background 
solution (pojPo, *i'0 7 ■Wo) by its Newtonian counterpart; 



the corresponding changes to 5p, 5p, 5v and 6ip are of 
post-2-Newtonian order and can be neglected. 

Equations ((5H|) . ((57T5)) and ((ETC)) together with the 
boundary condition (j5.10cp are precisely the standard 
Newtonian perturbation equations supplemented by the 
external acceleration a^xt- We have therefore shown that 
the leading order effect of the external gravitomagnetic 
tidal field Bij{t) on the star can be computed using New- 
tonian perturbation theory. The expression (|5.17p for 
the external acceleration a^xt can be rewritten using the 
formula (|5.12p for 5C as 



alxt = --e^jkBki{t)x'x^ -2e,jkVpki{t)xK (5.18) 

Using this expression together with the mode function 
(|5.7p and Vq = CI x x in the coupling integral (|2.5p writ- 
ten in inertial frame coordinates shows that the only non- 
vanishing driving terms occur for 1 = 2, and we obtain 



(^2m, lext) 



mniMR^I- 



15\/6 



2m ^imUt I^y2my-t 



Bij ^ 2Bpi^jpq^q 



15V6 



dt 



[e*'""*(3^2")*6„(0] . 

(5.19) 



Here are the dimensionless coupling parameters 

(|5.6p . and we have used H = (0, 0, fl). 

So far the analysis has been valid for a star placed in 
an arbitrary gravitomagnetic tidal field Bij{t). We now 
specialize to a star in a binary. We denote by z{t) and 
z{t) the coordinate location and velocity of the compan- 
ion star of mass M'; these quantities are the relative 
displacement and relative velocity since we are working 
in the center of mass frame of the star M. In Appendix 
|X]we show that 



Bij = &M'zi^iej)kiZkZi/r^. 
We parameterize the quasi-circular orbit as 

z{t) = r{t) [cos i/j cos (f>{t), sin (t){t), sin cos 0(i)]. 



(5.20) 



(5.21) 

Here (j){t) is the orbital phase of the binary, and ip is 
the inclination angle of the orbital angular momentum 
relative to the spin axis of the star. Inserting this pa- 
(5.17) rameterization into the formula (|5.20p for Bij and then 
into the expression (|5.19p for the overlap integral gives 
the following results for the I = 2,m = 1 and I = 2,m = 2 
r-modes 
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5 



(cj + f7) sin ( ^ ) sin ( ^ ) e*-^**) + {uj - 9.) cos ['^ ] cos [ ^ ] e 



(62,ae.t) = ^/^^MM'i?W2^'sinV' 



(2f] + c^) sin2 e'-^^*) - (2f7 - uo) cos^ (^^ 



i4,{t) 
(5.22a) 
(5.22b) 



where uj = (p. 

These driving terms are each a sum of two terms pro- 
portional to e""^*='=*'^(*) which oscillate at different fre- 
quencies. From the equation of motion (j2.4p . resonant 
driving will occur at uj{t) = luq if 



(5.23) 



Now the mode frequency uj^ is negative and by Eq. ([57 
satisfies < ^- Also the orbital angular velocity uj{t) 
is by convention positive. It follows that only the terms 
proportional to e™f^*-«0(*) can produce resonant driving. 
The other non-resonant terms oc ^ill produce 

corrections to the late time phase evolution that are at 
least a factor of e smaller than those due to the resonant 
terms, and can be safely dropped. Thus the resonant 
orbital frequency is 



LOq 



(5.24) 



directly use the results of section IIIII to obtain the reso- 
nant response of each mode. From Eq. (|3.22[) the quan- 
tities needed to characterize this resonant response are 
the resonant timescale iros [Eq- (|3.2p ] and the amplitude 
Aqj of the forcing term at resonance. From Eqs. (|5.26[) 
and (|5.25p we obtain for the values of these parameters 



ires,2i = [96A^^/'(2r!/3)^^/'75] 

^0,21 = 



-1/2 



(5.27a) 



— ^T^°HTr°H2 



'0 



ires,22 = [96A^ (4r!/3) 3/5 

^0,22 = 



-1/2 



167r /22 M'i? 417 . , 2 / V 

t: sm w cos — 

5 r2 3 ^ I 2 



(5.27b) 
(5.27c) 

(5.27d) 



D. Effect of current quadrupole on orbital motion 



Using the formula (|5.2p for the mode frequencies c^m = 
'-^im we find for the two modes we consider 



m 



m 



-n/3, 
-2n/3, 



Wo 



2f7/3, (5.25a) 
= 417/3. (5.25b) 



We next substitute the resonant terms from the overlap 
integrals (|5.22p into Eq. (|2.4p and use the formula (15. 
for b. This finally gives 



C21 + 



I nenhiM'R , 

2 '^\^ ~ ") 



5 

3^ 



X cos I — I cos I - 1 p^^-^m 



(5.26a) 



C22 + iW22C22 



^ bjyuj — 2\l) 



'^12 



X sin -0 cos ( ) 6 



"i- [ ^ \ „2int-i4>{t) 



(5.26b) 



Since r-modes do not generate mass multipole mo- 
ments to linear order in the Lagrangian fluid displace- 
ment^^, the Newtonian equation of motion for the star's 
center of mass worldline (and its companion's as well) 
remain unchanged when such modes are driven. The 
leading order correction to the equations of motion is 
a post-l-Newtonian tidal interaction term involving the 
current quadrupole moment induced in the star by the 
r-mode. The equations of motion including the effect of 
this current quadrupole were derived in Ref . [HBl • We spe- 
cialize Eq. (6.11) of Ref. [1^ to two bodies with masses 
Ml = M, M2 = M', positions Zi{t) and 22 (i), current 
quadrupoles S}j{t) — Sij{t) and Sfj{t) — 0, and with 
all other mass and current moments equal to zero. This 
gives 



M' , 
r3 Mt 



(5.28a) 
(5.28b) 



Note that for the aligned case = 0, only the Z = 2, to = 
1 mode is excited. No modes are resonantly excited in 
the anti-aligned case ^ = tt. 

The mode amplitude evolution equations (j5.26p are ex- 
actly of the form p. lip , with m' = 1. We can therefore 



This is correct only to leading order in the star's spin frequency. 
There are corrections to the leading order mode eigenfunction 
that scale as Q^R^/M. As discussed in the introduction, 
these correction terms can couple to the Newtonian tidal field and 
induce mass multipole moments to linear order in the Lagrangian 
fluid displacement |31| . 
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where z = Z2 — Zi, v = z, and the contribution a to the 
relative acceleration z due to the current quadrupole is 



M 



UpqSqr ^ 



^^ipq^j Sqr 



(5.29) 



Here the angular brackets denote taking the symmetric 
trace-free part of the corresponding tensor, thus 

^2 

Z<^iZj — ZiZjZf ~^^^ijZj' ^ Sij-Zj Sji-Zi^. (o.30) 

o 

In Sec. HI] we showed that only the tangential compo- 
nent of the relative acceleration is needed to compute 
the leading order phase lag parameter in the circular ap- 
proximation. That tangential component is given by 

1 



ai — — V ■ a 

Now for an axial fluid displacement ^{t,x), the induced 
current quadrupole moment to linear order in ^ is given 
by the following integral in inertial coordinates 



■'qr 



p{r) - ^k^kX<:qXr> + r'^^<q^r> 



qJ^r> 



(fx. (5.32) 



Using the mode expansion (|2.3p and the mode function 
(|5.3p we evaluate the integral (|5.32p for an Z = 2 r-mode. 
Only the term proportional to ^ survives, and we obtain 



Stt 2 



Sqr -— —\l -Ilm.MR 



(5.33) 

Substituting the current quadrupole (|5.33[) into the ac- 
celeration (|5.3ip and using the orbit parameterization 
(I5.2ip gives for the ra — 1 case 



ax. 
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x \ cos(3V/2)cos(V/2) [e"^-'^\c2i - inc2i)\ 



3in(3V'/2)sin(V'/2) [e 



(C21 



-l-C.C. 



(5.34) 



Next we use the fact that in the vicinity of the reso- 
nance we have Cj = ~iujjCj[l + 0(e)], from Sec. IIII 51 to 
eliminate the time derivative terms. We also drop the 
terms proportional to sin(3'(/'/2) sin('0/2) as they cannot 
give rise to secular contribution to the phase shift. This 
yields 

«^;2i = ^Tf /21 ^^""''Z M,R' cos(3^/2) cos(V'/2) 



C2l{t) 



c.c. 



(5.35) 



Similarly we find for the m = 2 mode 

[tT f2(w22 + 2fl) 



h 



MtR" sin?Acos^(V'/2) 



-le 



i(t>{t)-2iflt 



C22(t) 



(5.36) 



These tangential accelerations (|5.35l) and (|5.36p coincide 
with the prediction (j2.18bp of the Newtonian model of 
Sec. ini when we use the formulae (I5.27bp and (|5.27dp 
for the amplitudes at resonance ^o.j of the mode forcing 
terms. We can therefore use the results of this Newto- 
nian model. Substituting the amplitudes and resonance 
timescales (|5.27p into the formula (|4.12p for the phase 
shift A(j> and using Eqs. (|5.5p and (|5.25p finally yields 



A(l)2 



2/3 



5^ 
192 V3 

n'^/^R" 

M'M'^m}'^ ' 



^21 ' 



3?A 
"2" 



cos 



(5.37a) 



A0: 



'22 



5^ 
192 V3 



2/3 



I22 sin^ V cos^ 



M'M^M, 



1/3- 



(5.37b) 



E. 



Numerical values for specific neutron star 
models 



1. Barotropic stars 

In the case of barotropic stars (i.e. stars without buoy- 
ancy forces) the formula (|5.3p for the mode function only 
applies for |m| — I, in which case fi„i{r) oc r'. In this 
case, using the normalization condition (|5.4p . the formula 
(|5.6p for the parameter /22 evaluates to 



-^22 — 



1 



MR^ 



p{r) 



r^dr 



1/2 



This gives 



I22 — 



3 
28^ 



0.185 



(5.38) 



(5.39) 



for constant density stars, while for a F = 2 polytrope 
for which p{r) — Msin(7rr/i?)/(4i?^r) we obtain 



^22 — 



VI2O - 207r2 + TT^ 

2^^572 



0.128. 



(5.40) 



For other stable polytropes, one needs to solve the Lane- 
Emden equation numerically to obtain the number J22. 
For r = 4/3, 3/2 and 3 respectively, the number /22 turns 
out to be equal to 0.0437, 0.0812 and 0.155, a variation 
of about a factor of 5 over the range of all polytropic 
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indices. Evaluating the phase shift (j5.37bp for the equal 
mass case M — M' and at the value ip — tt/3 of the 
inclination angle which maximizes the mode driving gives 



A(/.22 = -0.05 Rtofs(ooM[.r^' 



for the constant density case and 



-0.03 RtofslooM 



10/3 
1.4 



(5.41) 



(5.42) 



for the r = 2 case. Here we have defined R = 10i?io km, 
M = IAMiaMq, and n/{2n) = 100/siooHz. As dis- 
cussed in the introduction, these phase shifts may be 
marginally detectable by gravitational wave interferome- 
ters for nearby inspirals. 

Consider next the case m ^ I. As discussed above 
there are no purely axial modes of this type in barotropic 
stars. There are modes analogous to the r-modes which 
have both axial and polar pieces; these are the inertial or 
hybrid modes [6l|. These modes are somewhat difficult 
to compute in barotropic stars, but are easier to com- 
pute in stars with buoyancy where they are purely axial. 
Therefore for the m ^ I case we will switch to using a 
slightly more realistic neutron star model which includes 
buoyancy. 



2. Stars with buoyancy 



In stars with buoyancy, the formula (j5.3p for the mode 
functions is valid for all m and /. However in this case 
the radial mode function fim(i') is not a simple power 
law, but instead must be obtained from solving numer- 
ically a complicated Sturm-Liouville problem 36] . The 
stellar model we use consists of a F = 2 polytrope for 
the background star, with the perturbations character- 
ized by a constant adiabatic index Fi ^ F chosen so 
that the Brunt- Vaisala frequency evaluated at a location 
halfway between the center of the star and its surface 
is 100 Hz. See Refs. HHl, [fil for detailed models of 
the buoyancy force and corresponding ij-modes of which 
our chosen model is representative. Our results below 
are fairly robust with respect to the choice of buoyancy 
model, varying by no more than ^ 20% for a ~ 50% 
change in the Brunt- Vaisala frequency evaluated at a lo- 
cation halfway between the center of the star and its 
surface. 

We compute the radial mode function for modes 
with different numbers n of radial nodes. The net phase 
shift is obtained by summing over all these modes, so we 
use an effective value of /|i obtained by summing over 
n. We find that it is sufficient to sum over n = 0, 1 and 
2 in order to obtain two digits of accuracy in , and we 
obtain 



Jo 



0.016. 



(5.43) 



We substitute this into Eq. (|5.37ap and evaluate the 
phase shift for the equal mass case, and at the value 



■0 = of the inclination angle which maximizes the mode 
driving. This gives 



A02 



-0.03 RUlimM,l°'^ 



(5.44) 



VI. CONCLUSIONS 



In this paper we have studied in detail the gravitomag- 
netic resonant tidal excitation of r-modes in coalescing 
neutron star binaries. We first analyzed the simpler case 
of Newtonian resonant tidal driving of normal modes. 
We showed, by integrating the equations of motion of 
the mode-orbit system, that the effect of the resonance 
to leading order is an instantaneous shift in the orbital 
frequency of the binary at resonance. The effect of the 
modes on the orbital motion away from resonance is a 
small correction to this leading order effect. We then 
showed that for the case of gravitomagnetic tidal driv- 
ing of r-modes, the equations of motion for the mode- 
orbit system can be manipulated into a form analogous 
to the Newtonian case. We then made use of the New- 
tonian results to compute the instantaneous change in 
frequency due to the driving of the r-modes. This shift 
in frequency is negative, corresponding to energy being 
removed from the star and added to the orbital mo- 
tion. This sign of energy transfer occurs only for modes 
which satisfy the Chandrasekhar-Friedman-Schutz insta- 
bility criterion. We also argued that the resonances are 
at best marginally detectable with LIGO. 
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APPENDIX A: COMPUTATION OF THE 
GRAVITOMAGNETIC TIDAL TENSOR 

In this appendix we compute the gravitomagnetic tidal 
tensor Bij{t) using two different methods, first using the 
global inertial frame, and second directly in the star's 
local asymptotic rest frame. 

We denote by (?,£*) the global harmonic center-of- 
mass frame. In this frame the Newtonian potential and 
gravitomagnetic potential are 



$(i,i) 
C(i,i) 



M 



M' 



\A-z,{t)\ 
AMzi{i) 



\x-zi{t)\ \X-Z2{t)\ 



(Ala) 
(Alb) 
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The post-Newtonian scalar potential -0 will not be needed 
in what follows. We now transform to a body-adapted 
coordinate system {t,x^) of the star of mass M. We 
use the general coordinate transformation from one har- 
monic, conformally Cartesian coordinate system to an- 
other given by Eq. (2.17) of Ref. [HI] with {t, x^) replaced 
by {t,x'^) and {t,x^) replaced by (t, a;*). That coordinate 
transformation is parameterized by a number of free func- 
tions: a function ac (t) governing the normalization of the 
time coordinate at Newtonian order; a function (P'{t) gov- 
erning the translational freedom at Newtonian order (de- 
noted in Ref. [56] ) ; a function Rk (i) governing the an- 
gular velocity of the coordinate system; a function hl{t) 
governing the translational freedom at post-Newtonian 
order, and a free harmonic function /3c{t,x''). The first 
four of these functions are fixed by the requirements listed 
at the start of Sec. IV Al above. The fifth function (3c pa- 
rameterizes the remaining gauge freedom, under which 
the tidal tensors are invariant. In fact the gravitomag- 
netic tidal tensors do not depend on hi, so we can use 
a coordinate transformation with h^. — Pc — 0. The re- 
quired values of the other functions are given by 



d{i) 
R{t) 



Mi), 

X Ccxt[i,2i(?)] + 2ii(t) X V$ext[i,2i(i)] 



(A2a) 
(A2b) 



+ ^Zi{i) X ii(t). 



(A2c) 



where $oxt and (^cxt denote the second terms on the right 
hand sides of Eqs. (|Al|l . 

Using this coordinate transformation we can compute 
the transformed gravitomagnetic potential C{t, x) in the 
body frame (t, a;*). We can then extract the gravitomag- 
netic tidal moment Bij(t) = —Hij(t)/2 by comparing 
with the general multipolar expansion of C given in Eq. 
(3.5c) of Ref. [1^. That extraction method is equivalent 
to throwing away the pieces of C that diverge at the ori- 
gin a; = and evaluating Bij = —B^i^j^ at a; = 0, where 
B = V X is the gravitomagnetic field. After eliminating 
the acceleration terms using the Newtonian equations of 
motion the result is the expression (|5.20p . 

An alternative, simpler method of computation is to 
solve the post-l-Newtonian field equations directly in 
the body frame. The only complication here is that the 
boundary condition on the potentials $ and ^ as r ^ oo 
are no longer $ ^ 0, ^ ^ 0, since the body frame is ac- 
celerated and rotating with respect to the global inertial 
frame. In the body frame the two stars are located at 
a; = and at x = z{t) = Z2{t) — Zi{t). The solutions to 
the field equations are thus 



^t,x) = ---- 



M' 



C{t,x) 



x\ \x — z{t)\ 
AM'z{t) 



\x~z{t)\ 



+ $bc(i,a;), (A3a) 
Chc{t,x). (A3b) 



Here $bc and ^bc are solutions of the Laplace equation 
that are determined by the large-r boundary conditions. 
They are determined by the functions ac, (i% h]. and R'' 
discussed above. The key point is that these terms do 
not contribute to the gravitomagnetic tidal tensors; this 
can be seen from the explicit form of these terms given in 
Eq. (3.38c) of Ref. [sH]. Therefore we can drop the terms 
^bc and Cbc and extract the gravitomagnetic tidal tensor 
Bij{t) using the method discussed above; the result is 
again given by Eq. (|5.20p . 



APPENDIX B: ALTERNATIVE DERIVATION OF 
TIME DELAY PARAMETER USING 
CONSERVATION OF ENERGY 

In this appendix we give an alternative derivation of 
the time delay parameter At using conservation of en- 
ergy, for the Newtonian model of Sees . HIl Hill and HVl and 
also for the r-mode driving of Sec. |Vl The energy con- 
servation method is simpler to use than the method of 
integration of the orbital equations of motion used in the 
body of the paper. However, it only gives information 
about the time delay parameter Ai, and says nothing 
about the phase shift parameter A^. In order to de- 
rive the relation (14. lip between these two parameters, an 
analysis of the orbital equations of motion is necessary. 



1. Newtonian mode driving 

Consider the evolution of the binary from some initial 
frequency Ui below resonance to some final resonance Uf 
after resonance. We denote by 



U!{uj) — UJ-p-p{uj) + 5lli{uj) 



(Bl) 



the time derivative of the orbital frequency w as a func- 
tion of L0\ this consists of a sum of the inspiral rate LO-pp[ijj) 
for point particles [Eq. (|1.5p above] together with a per- 
turbation 56j{u!) due to the mode coupling. The total 
time taken for the binary to evolve from uji to w/ can be 
written as 



duj 

uj{uj) 
duj 



UJr. 



>(^) 



6lu{uj) 

^pp(w)^ 



duj. (B2) 



Here on the second line we have linearized in 6lu. The 
first term in Eq. (|B2p is the time it would take for point 
particles, and by comparing with Eq. (|4.10p we can iden- 
tify the second term with the negative of the time delay 
parameter At: 



At = 



6u!{lu) 



dui. 



(B3) 



Consider next the total energy radiated into gravita- 
tional waves between oji and LOf. This radiated energy 
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can be written as 

E,,4cu,^ujf) = r ^^duj (B4) 

(B5) 

where Eg„{uj) is the gravitational wave luminosity and 
as before we have linearized in A key point now 

is that the gravitational wave luminosity Egvj{^) is to a 
good approximation unaffected by the mode coupling, so 
that it is the same function for the point-particle inspiral 
and for the true inspiral. "'^^ Therefore the first term in 
Eq. (|B5p is just the energy radiated between oji and LOf 
for a point-particle inspiral, namely 

Here Eo^\y{u)) is the orbital energy of the binary at fre- 
quency OJ. 

Consider now the evaluation of the second term in Eq. 
(IBSp . The function 5uj(ijj) is peaked around the resonance 
Lj — ljJq with width Aw ~ tj(wo)^rcs ~ cwq. The fractional 
variation in the luminosity Eg^(uj) over this frequency 
band is of order ~ e, so to a good approximation we 
can replace £'gw(w) with Eg^^LOo) and pull it outside the 
integral. "'^'^ The remaining integral is just the formula 
(|B3p for the time delay parameter At, and we obtain 

Erad{i^i ^ ^f) = -^orb (^^i ) ~ -^orb (l^/ ) 

-Eg„{LOo)At[l + 0{e)]. (B6) 
Finally global conservation of energy gives 

^^orb(Wi) -|- -E'star(Wi) = £'rad (l^i ^ ^f) + Eoih{^f) 

+E,t,r{iOf), (B7) 

where i?star is the total energy of the star.^^ Defining 

AE'star — Estar{iOf) — Egtari^i) (B8) 



and combining this with Eq. (|B6[) yields 

Ai^star = 4wK)Ai [1 + 0(e)] . (B9) 

Note that a simple interpretation of the formula (|B6[) 
is that the resonance can be treated as an instantaneous 
change of frequency of magnitude Auj = Wpp(a'o)Ai and 
corresponding instantaneous change in orbital energy of 
Egvj{LOo)At. However this interpretation is not com- 
pletely valid since the resonance is not instantaneous: 
the energy radiated during resonance is larger than the 
energy Emodc absorbed by the mode by a factor ~ 1/e. 
The formula (jB6|) must therefore be derived by integrat- 
ing over the resonance as we have done above. 

The change in the energy of the star AEgtai consists 
of the rotating-frame energy in the mode i^modo, the re- 
maining portion of the energy in the mode due to its 
angular momentum, and the change in the rotational ki- 
netic energy of the star. The total change can be com- 
puted from a prescription detailed in Appendix K of Ref. 
[irt . which gives 

AE,t.r = ^^[uj.+mjQ]. (BIO) 

LUj 

Here ^mod c is given by Smodo = bjUJj\cj\'^ = bjUjjBf, cf. 
Eq. (|2.16p above. Combining this with Eq. (jB9p gives 
a formula for At which agrees with our earlier formula 
(|4.9p . Equation (|B10p is computed as follows. Equation 
(K17) of [ZT] gives the formula for the physical energy 
due to the fluid perturbation. The first term of (K17) is 
the rotating-frame energy of the mode and the two other 
terms are angular momentum contributions to the total 
energy of the perturbation. The rotating-frame energy is 
given by (K21) and the angular momentum contribution 
is computed using the binary's orbital equations of mo- 
tion including mode coupling, along with conservation of 
total angular momentum following the work of Ref. 

We conclude that using energy conservation gives a 
simple method to evaluate the time delay parameter At, 
without computing the evolution of the orbital variables. 
However, energy conservation does not give any direct in- 
formation about the phase shift parameter Acf); for that 



^ This approximation is valid for computing the resonance phase 
shift. It is not vaUd, however, for computing the much smaller 
orbital phase correction in the adiabatic regime before resonance. 
In that regime one must include the gravitational radiation from 
the mode, which is phase coherent with that from the orbit [l9l 

m. 

^ Note that this approximation can only be applied to the second 
term in Eq. IIBSI l and not to the entire integral ||B4| I: the corre- 
sponding errors in energy would in that case be larger than the 
energy we are computing here. 

* We neglect here the tidal interaction energy corresponding to the 
last two terms in the Hamiltonian 1 12.161 1. Although this energy 
is comparable to the mode energy during resonance, at late times 
after resonance it is smaller than the mode energy by a factor 
~ e. 



Note that the splitting of the total energy of a perturbed star 
into mode and rotational contributions depends on the choice of 
a background or reference uniformly-rotating star. In the for- 
malism of Ref. 47], such a choice is determined by the use of 
Lagrangian displacement as the fundamental variable. Other 
choices such as Eulerian velocity perturbation yield different re- 
sults for this splitting. The same issue applies to angular mo- 
mentum, and explains why the angular momentum of r-modes 
was found to be zero in Ref. |64| , in disagreement with the results 
of Ref. [13. 

Over the resonance timescale, the amount of angular momen- 
tum radiated away by gravitational waves is negligible for this 
computation. 



22 



one must analyze the resonance itself. As we have dis- 
cussed, the result of this analysis is the relation (|4.1ip be- 
tween and At, and using this relation together with 
Eq. (|B9p we obtain a formula for the phase shift in terms 
of the mode energy: 



A(j) = 



mjUJo\ Wo-E'modo 



(Bll) 



In Ref. [3l| energy conservation was used to compute 
the phase shift, although there At/) was misinterpreted as 
the asymptotic value of the phase perturbation after the 
resonance. The fact that the phase perturbation grows at 
late times [cf. Eq. p.9p above] was noted in Refs. j29l . [6^ . 



2. Gravitomagnetic mode driving 

The argument of the previous subsection can be carried 
over to post-l-Newtonian gravity with a few minor mod- 
ifications. First, the total conserved post-l-Newtonian 
energy E can be split up into two pieces, an "orbital" 
piece Eoih and a remaining piece, as follows. On a given 
constant time slice we can compute the post-l-Newtonian 
positions and velocities of the stars' center-of-mass world- 
lines, and also the post-l-Newtonian masses of the bodies 
[sgI \5T\ . The post-l-Newtonian masses will evolve with 
time, so we choose the values of these masses at t — > —oo 
when the stars are far apart. We then compute the to- 
tal post-l-Newtonian energy for a system of point par- 
ticles with the same positions, velocities and (constant) 
masses; this defines the orbital energy E'orb of the sys- 
tem. The remaining piece of the total energy consists of 
the mode excitation energy -Emodo, together with an in- 
teraction energy i^int oc StjBij which we can neglect for 
the same reason as in the Newtonian case.^^ With these 
definitions the arguments of the last subsection are valid 



for gravitomagnetic driving. 

At late times it is smaller than -Emodc by a factor of e. 



In the Newtonian case it was most convenient to eval- 
uate the energy i^mode transferred into the stellar modes 
by using the explicit formula (jBlOp for the mode energy. 
In the post-Newtonian case it is more convenient to use a 
different technique. A formula for total tidal work done 
on the star by the gravitomagnetic driving can be derived 
by surface integral techniques [67| : the result is 



dE, 



^ = -|6.,(i)^.,(i). 



(B12) 



See Refs. [H, [6^ [t^, [zH for further discussion of the 
derivation method and of similar tidal work formulae. 



We now use the explicit formula (|5.20p for the grav- 
itomagnetic tidal moment Bij, and compare the result- 
ing expression with the formula ()5.3ip for the tangential 
acceleration due to the current quadrupole moment, 
specializing to a circular orbit. The result is 



dt 



(B13) 



Using the formula (|B9p for At and integrating over the 
resonance gives 



At: 



£'gw(wo) 



dtaM. 



(B14) 



This is the same result as was obtained in Sees. IIV Bl and 
IVDI for the time delay parameter; those sections also in- 
volved integrating with respect to time, and it can 
be checked that the prefactors match. Thus energy con- 
servation gives the same result as the method of direct 
integration of the equations of motion in the gravitomag- 
netic case. 
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